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Abstract 

We propose the model, which allows us to approximate fractional Levy noise 
and fractional Levy motion. Our model is based (i) on the Gnedenko limit 
theorem for an attraction basin of stable probability law, and (ii) on regarding 
fractional noise as the result of fractional integration /differentiation of a white 
Levy noise. We investigate self - affine properties of the approximation and 
conclude that it is suitable for modeling persistent Levy motion with the Levy 
index between 1 and 2. 
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I. INTRODUCTION. 



By Levy motions, or Levy processes, we designate a class of random functions, which 
are a natural generalization of the Brownian motion, and whose increments are stationary, 
statistically self-affine and stably distributed in the sense of P. Levy JE]. Two important 
subclasses are (i) the stable processes, or the ordinary Levy motions, which generalize the 
ordinary Brownian motion, or the Wiener process, and whose increments are independent, 
and (ii) the fractional Levy motions, which generalize the fractional Brownian motions and 
have an infinite span of interdependence. 

The theory of processes with independent increments was developed beginning from the 
Bachelier's paper || concerning Brownian motion. However, the rigorous construction of 
this process and studies of properties of its trajectories were undertaken by Wiener [Q. 
The modern presentation of the general theory of processes with independent increments is 
contained in ||||. The theory of stable processes with independent increments has begun 
its history from the already cited work |I| and, later on, was developed by other prominent 
mathematicians. In particular, the properties of extrema of stable symmetric processes were 
studied in Ref. ||. The local properties of stable processes and geometric properties of their 
trajectories were considered in Ref. 0. 

The basis of the general theory of random processes with stationary self-affine increments 
was laid by Kolmogorov ||. Fractional Brownian motions were introduced by Mandelbrot 
and van Ness as a (relatively) simple family of random functions "that could be in some way 
be expected to be "typical" of what happens in the absence of asymptotic independence " 
||. In this paper the possibility was also pointed out of constructing "fractional Levy-stable 
random functions" in a way analogous to that of constructing fractional Brownian motion. 
Such stable generalizations of fractional Brownian motion were introduced in p!0|JTT| . 

The Levy random processes play an important role in different areas of application for 
at least two reasons. 

The first one is that the Levy motion can be considered as a generalization of the Brow- 
nian motion. Indeed, the mathematical foundation of the generalization are remarkable 
properties of stable probability laws. From the limit theorems point of view, the stable dis- 
tributions are a generalization of widely used Gaussian distribution. Namely, stable distri- 
butions are the limit ones for the distributions of (properly normalized) sums of independent 
identically distributed (i.i.d.) random variables fl2| . Therefore, these distributions (like the 
Gaussian one) occur, when the evolution of a physical system or the result of an experi- 
ment are determined by the sum of a large number of identical independent random factors. 
An important distinction of stable probability densities is the power law tails decreasing as 
, x — > oo , a is the Levy index, < a < 2. Hence, the distribution moments of the 
order a diverge. In particular, stably distributed variables possess a non-finite variance. 

The second reason for ubiquity of the Levy motions is their remarkable property of scale - 
invariance. From this point of view the Levy motions (like the Brownian ones) belong to the 
so - called fractal random processes. Indeed, the objects in nature rarely exhibit exact self - 
similarity ( like the von Koch curve), or self - affinity. On the contrary, these properties have 
to be understood in a probabilistic sense fll3Ul4j|. The random fractals are believed to be 
widely spread in nature. A coastline is a simple example of statistically self - similar object, 
whereas the trace of the Brownian motion is statistically self - affine. Several numerical 
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algorithms were developed in order to simulate fractional Brownian motion p5|Jl^J16|J17 . 
They allow one to model many highly irregular natural objects, which can be viewed as 
random fractals ||14||. The traces of the Levy motions are also statistically self - affine, 



therefore, one may expect that they are also suited for modeling and studies of natural 
random fractals. 

The stable distributions and the Levy, or Levy - like, random processes are widely used 
in different areas, where the phenomena possessing scale invariance (in a probabilistic sense) 
are viewed or, at least, can be suspected, e.g., in economy fT8Lp0f, biology and physiology 
PH , turbulence |Z| and chaotic dynamics P3| , solid state physics pi| , plasma physics p5| , 
geophysics |26| etc. In this respect, the models are needed, which allow one to simulate 
fractional Levy motion (fLm) with the prescribed quantitative statistical properties and to 
test and improve methods aimed at analysis and interpretation of experimental data. As 
far as the authors know, the properties of fLm are, at least, rarely discussed in physics 
literature. Therefore, we believe, that the models constructed "at the physical level of 
strictness" may be useful for developing effective methods for experimental data processing 
and for the purposes of numerical modeling. When studying this problem, we find that 
the different levels of complexity are required for the models aimed to approximate fLm, 
which differ by their Levy index a and by their " long- memory" behavior, that is, persistent 
or anti-persistent one. We discuss this item below and restrict ourselves with the simplest 
model, which serves as a good approximation to the persistent fLm with the Levy index 
lying between 1 and 2. We propose a simple approximation to persistent fLm and study its 
scaling properties. 



II. "r H LAWS" FOR THE FRACTIONAL LEVY MOTION. 

Our aim is to make an approximation to the process denoted by L a ^(t), whose incre- 
ments are stably distributed with the Levy index a, < a < 2, and possess the properties of 
stationarity and self - affinity By analogy with the Definition 3.2 of Ref. M, the increments 
of a random function L a ^(i) will be said to be self - affine with parameter H if for any 
k > and any t 

L a , H (t + KT) - L a>H (t ) = K H (L ajH (t + T) - L ajH (t )) (1) 

The particular cases are: 

(i) ordinary Brownian motion, L a ji(t) = B(t), which has H — 1/2; 

(ii) fractional Brownian motion, L a> H{t) = Bn{t), which has < H < 1; 

(iii) ordinary Levy motion, L aj #(t) = L a (t), which has H = 1/a. 

Since the Levy-stable distribution with the Levy index a possesses the moments of the 
order q < a, the "t h laws" (by terminology of |§) for the structure functions of the fractional 
Levy motion can be stated as follows: for any < q < a 

(\L a , H {t + T)-L a , H {t)\ q ) 1,q ocr" (2) 

where the proportionality coefficient is a function of q and the parameters of stable distri- 
bution for the increments, but not of H . The range of H can be determined as follows. By 
Minkowski's inequality, for any positive r 1; r 2 > and 1 < q < a 
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or, using Eq.(2), 



(lW(i + ri+r 2 )-L tt)fl (t)f) 1/? < 
\L a , H (t + n + r 2 ) - L a , H (t + T!)| 9 ) 1/9 + 

+ (iw(i+n)-wwn 1/? (3) 

(T! + r 3 )" < Tf + rf (4) 



which implies H < 1. Further, we require g - th order mean continuity of L ai #(i), that is [|J, 
for any t and g < a 

^{\L ehH {t + M)-L ehH {t)\ t )=0 (5) 

It implies that H > 0. 

From Eq. (1) the "r^ law" follows also for the range of the fLm with the Levy index 
a > 1, 

R(r) = sup (Lo,ff(s) - L QtH {0)) - inf (L Q , H (s) - L a , H (0)) (6) 

0<S<T 0<S<T 

(i2(r)> oc r H (7) 

We make some remarks concerning distinction between the motions with the Levy index 
lying between and 1, and those with the index lying between 1 and 2. Comparatively to 
the latter case, the former one is more involved for the analysis and numerical simulation. 
The reasons are as follows: 

(i) since q must be less than a, the Minkowski's inequality does not hold for a < 1, hence 
we can not use it for establishing the upper boundary for H; 

(ii) the statistical mean of the range is infinite for a < 1; 

(iii) with a decreasing, the increments of fLm grow rapidly due to the power - law tails 
of stable distributions. Very large increments are occurring more and more frequently, and 
the moments of the increments strongly fluctuate from realization to realization. Thus, 
the reliability of the results of simulation decreases. The same conclusion can be done 
when treating experimental data, in which the Levy statistics with a < 1 can, at least, be 
suspected. 

The three peculiarities mentioned above require a modification of the methods, which 
we use below for making an approximation to fLm and for studying its properties. For this 
reason, we restrict ourselves by the case 1 < a < 2, however, having in mind to study more 
complicated case < a < 1 in future. Fortunately for us, there are evidences from different 



areas of application, see, e.g., |L8| [Q, that, as a rule, the Levy processes with a > 1 do 
occur. 
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III. THE MODEL. 



The process of constructing approximation to fLm can be divided into 3 steps. 

Step 1. At the first step we generate random sequence of i.i.d. random variables possess- 
ing stable probability law. These variables play the role of increments of the ordinary Levy 
motion having the Levy index a, < a < 2. The value a = 2 corresponds to the ordinary 
Brownian motion, hence in this case the sequence of independent increments is generated 
with the use of a standard Gaussian generator. Since in [[OJ the sequence generated at 
a = 2 is called " approximate discrete-time white Gaussian noise" one may call the sequence 
generated at < a < 2 " approximate discrete-time white Levy noise" . 

We restrict ourselves by symmetric stable laws with the probability density p a ,D(x) an d 
the characteristic function 

Pa,D{k) = (exp(zfcc)) = exp(— D \k\ a ) . (8) 

Here a is the Levy index, and D is a positive parameter. At a = 1 and 2 one has the Cauchy 
and the Gaussian probability laws, respectively. In other cases the symmetric stable laws 
are not expressed in terms of elementary functions. At < a < 2 they have power law 
asymptotic tails P"2"f, 



T(l + a)sin(7ra/2) 

Pa,D[X) OC D ■ TjTCJ , X -> ±00 . (9) 

Among the methods of random sequence generation with the given probability law F(x) the 
method of inversion seems most simple and effective [|27||. However, it is well-known that 



its validity is limited by the laws possessing analytic expressions for hence, the direct 
application of the method of inversion to the stable law is not expedient. In this connection, 
we exploit an important property of stable distributions. Namely, such distributions are 



limiting for those of properly normalized sums of i.i.d. random variables [12|]. To be more 
concrete, we generate the needed random sequence in two steps. At the first one we generate 
an "auxiliary" sequence of i.i.d. random variables {£/}, whose distribution density F (x) 
possesses asymptotics having the same power law dependence as the stable density with the 
Levy index a has, see Eq.(9). However, contrary to the stable law, the function F(x) is 
chosen as simple as possible in order to get analytic form of For example, 

2 1 + \x\ 



Fix) = 1 r X > . 

v ' 2 (l + x a ) 



Then, the normalized sum is estimated, 



i m 

X =^Zfi • (ID 



where 
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\ 1/a 
7T \ 



" V2r(a)sin(7ra/2) J 

According to the Gnedenko theorem on the normal attraction basin to the stable law W2 



the distribution of the sum (10) converges to the stable law with the characteristic function 
(8) and D = 1. It is reasonable to generate the stable sequences with D = 1, doing the 
rescaling after that, if necessary. Repeating N times the above procedure, we get a sequence 
of i.i.d. random variables X(t),t = 1,2,..., N. This is an approximation to a discrete - time 
white Levy noise. 

Step 2. At the second step we convert approximate white Levy noise X(t) into the 
approximate fractional Levy noise X u (t). For this purpose one can use the procedure, which 



is called fractional integration/differentiation of a white noise |L5|. We remind the relation 
between the Fourier transforms of the function X(t) and of its fractional integral/derivative 
X u (t) of the u-th order: 

where v is positive in case of fractional integration and negative in case of fractional dif- 
ferentiation (to be more accurate, we say about left-side Riemann - Liouville fractional 
integral/derivative at the infinite axis |28|). Fractional integration of a white noise leads to 



amplification of low-frequency spectral components of the noise and thus, to the persistent 
process, whereas fractional differentiation leads to amplification of high-frequency spectral 
components and thus, to the anti-persistent process. At a qualitative level the difference 
between the persistent and anti-persistent behaviors can be formulated as follows: in the per- 
sistent random process the available tendency is supported, whereas in the anti-persistent 
process the opposite tendency prevails (T3| |2"5[| . The range of v will be discussed below. 
The approximate discrete - time fractional Levy noise X u (t) plays the role of the sequence 
of increments of the approximation to fLm. 

Step 3. With using approximate discrete - time fractional Levy noise X u (t) the approx- 
imation to the fractional Levy motion is defined by 

L a , v {t) = Y^X v {r) . (12) 

T=l 

Let us discuss the restrictions on the possible values of v. They are dictated by the 
restrictions on the range of H for the fLm, which we intend to approximate. 

At first, we note that the approximation to the white Levy noise possesses the property 
of self - affinity, 

X(Kt) = K l ' a X(t) . (13) 
Here k is a positive integer. It follows from Eq. (13) that 

X (^U ^ X M , (14) 
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and thus, according to our way of modeling, 



X v {Kt) = K U+1/a Xv{t) ■ (15) 

This is the property of self - affinity of the increments of the approximation to fLm. By 
comparing it with Eq.(l), we conclude that the equality 

H = v + l/a (16) 

must hold. Taking into account the restrictions on H, we conclude that for 1 < a < 2 

- 1/a < v < 1 - 1/a . (17) 

The particular cases of the proposed approximation £1X6 cLS follows: 

(i) a = 2, v = 0. In this case we get approximation to the Gaussian process with inde- 
pendent increments, that is, to the ordinary Brownian motion, which has the parameter H 
equal 1/2; 

(ii) a = 2, —1/2 < v < 1/2. In this case we get approximation to the fractional 
Brownian motion, which has the parameter H lying between and 1. We studied the 



properties of this approximation in [30 



(iii) < a < 2, v = 0. In this case we get approximation to the ordinary Levy 
motion, which has the parameter H equal 1/a. The step 2 is omitted when getting this 
approximation. Its properties were studied in JJT]. The three different models of the ordinary 



Levy motions were proposed in p2[ . They can be interpreted as a "difference scheme" to 
approximate the evolution equation for the density distribution of the ordinary Levy motion. 

In Fig.l the admissible range of the index v is shown for 1 < a < 2 on the (a, v)- 
plane. This range is bounded by the curve H = u + l/a = lm the top and by the curve 
H = in the bottom. The right vertical boundary indicated by thick line corresponds to 
the fractional Brownian motion, a = 2, — 1/2 < v < 1/2. The horizontal thick line v = 
corresponds to the ordinary Levy motion. This line divides two regions: upper one, v > 0, 
in which the Levy motion is persistent, and the bottom one, v < 0, in which the motion is 
anti-persistent. Dotted lines a, b relate to Figs. 4, 5 and will be explained below. 

Before proceeding with the numerical results illustrating the self-affinity properties of 
our approximation, we discuss the peculiarities of simulating persistent and anti-persistent 
motions. The problem appearing can be explained by taking fractional Brownian motion 
as an example |30| . At first we consider the persistent case. Fractional integration of a 
white Gaussian noise leads to the fractional Gaussian noise with spectral density decreasing 
as frequency increases. Thus, at numerical simulation the high-frequency inaccuracies are 
small and the low-frequency ones play the main role. Now we turn to the anti-persistent 
case. Fractional differentiation of a white Gaussian noise leads to the fractional Gaussian 
noise with spectral density increasing as frequency increases. At numerical simulation the 
high-frequency inaccuracies are prevailing. We have demonstrated, when studying fractional 
Brownian motion approximation based on fractional integration/differentiation of a white 
noise [30], that the high - frequency inaccuracies are more essential than the low - frequency 
ones, that is, the anti - persistent case is modeled with less accuracy. Moreover, it is known 
that the other algorithms for simulating fractional Brownian motion, namely those proposed 
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nj and in [14 as well, have the same drawback. Thus, one may expect that the anti- 



persistent fLm is also worse modeled. This is indeed so, as we convinced during numerical 
simulation. When simulating fractional Brownian motion, the problem is overcome in part 
by passing from "Type -2 Approximation" to a more sophisticated "Type 1 Approximation" 
having a finer grid jn|. A similar improvement can be performed when doing fractional 
differentiation of a white Levy noise. However, such a procedure requires more analysis, 
which have to be the subject of a separate paper. That is why below we present the results 
for the persistent case only. We also note that the persistent behavior is prevalent in nature 
|T5[ | |29| |, thus, we may hope that approximation to the persistent fractional Levy motion 
is more needed for applications. Of course, this does not imply that there is no need for 
studying anti - persistent fLm. On the contrary, we have in mind that the complex systems 
with the feedback require the development of approximations well-suited for modeling anti 
- persistent behavior. As an example we quote long - range anti - correlations and non - 
Gaussian behavior of the heartbeat of the healthy subjects 



IV. NUMERICAL RESULTS. 

The results of numerical simulation and analysis are shown in Figs. 2-5. 

In the top of Fig. 2 the probability densities p(x) for the members of the sequence X(t) 
are depicted by black points for (a) a = 1.2, and (b) a = 1.7. We use m = 30 terms in the 
sum (10). The functions p a ,i{x) obtained by the inverse Fourier transform, see Eq.(8), are 
shown by solid lines. In the bottom of Fig. 2 the black points depict asymptotics of the 
same probability densities in log - log scale. The solid lines show the asymptotics given by 
Eq. (9). It is seen that the Levy index can be estimated with the use of those values of 
X(t), which lie outside the peak located around x = 0. For the comparison the asymptotics 
of the Gaussian distribution with zero mean and unit variance are also shown. The log - log 
scale allows one to demonstrate smallness of probability of extreme Gaussian events. The 
examples presented demonstrate a good agreement between the probability densities for the 
sequences X(t) obtained with the use of the numerical algorithm proposed and the densities 
of the stable laws. 

We would like to note that the simplicity is a certain merit of the proposed approximation 
to a white Levy noise. The approximation is entirely based on classical formulation of 
one of the limit theorems and can be easily generalized for the case of asymmetric stable 
distributions. It also allows one, after some modifications, to speed up the convergence to 
the stable law. These problems, however, ought to be the subject of a separate paper. We 
also note, that two schemes were proposed recently, which use the combinations of random 
number generators |34j and the family of chaotic dynamical systems with broad probability 



distributions [pq| , respectively. 

In Fig. 3 typical samples of approximation to discrete - time Levy noises with the Levy 
index a = 1.5 are depicted for a white noise, v = (at the left), and for a fractional 
one, v = 0.3 (at the right). In both cases large frequent "outliers" are clearly seen, which 
are totally absent in case of the Gaussian noises. An interesting feature of the top figures 
is that they look very similar. Indeed, only careful comparison can reveal the differences 
between them. However, the corresponding ordinary Levy motion shown below at the left, 
and fractional persistent Levy motion shown below at the right are strikingly different. It 
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is so, because the fractional noise possesses such correlations, that the values of the noise 
at different instants " acts coherently" . It means that in case of persistent motion, not only 
"jumps", or "Levy flights" (which are clearly visible on both trajectories) lead to a large 
departure of the trajectory from the x - axis, but also "coordinated action" of the noise at 
different instants results in such a departure. In terms of anomalous diffusion we may argue 
that, contrary to the case of the ordinary Levy motion, when an anomalous diffusion rate 
is determined almost solely by the "Levy flights", in the persistent Levy motion it is also 
determined by "small steps". This is a qualitative corollary of the procedure of fractional 
integration of a white noise. 

In Fig. 4 we illustrate the properties of structure functions. We study r- dependence of 
the structure functions of our approximation, 

(\L^(t + T)-L a ^t)\ (1 ) 1/q ocr s , (18) 

where q < a. Since for the fLm the "t 11 law" is fulfilled, see Eq.(l), it is expected, according 
to our way of modeling, that s is close to v + 1/a. We demonstrate the results for q = 1/4. 
However, we verified that for any q less than a and not very close to it the results are just 
the same. We present the results for the approximation to the persistent Levy motions, 
whose parameters a, v vary along dotted lines a, b shown in Fig. 1. 

In Fig. 4a the results for the vertical line a are presented, a = 1.7. The exponent s versus 
v, see Eq.(18), is depicted by black squares. The expected relation s = v + 1/a is shown 
by solid line. We see that numerical results are well fitted by the expected line, but the 
discrepancy appears when v reaches its upper boundary, that is, for the strongly persistent 



case. The same effect appears when simulating fractional Brownian motion, a = 2 |30 
Since the Fourier transform of a white noise is divided by u u , v > 0, when getting fractional 
noise in the persistent case (see Step 2 of the approximation), one may suppose that the 
appeared discrepancy is due to the growth of low - frequency inaccuracies with v increasing. 
However, this is not a single reason: indeed, we convinced that making longer sample path 
does not diminish discrepancy essentially. 

In Fig. 4b the results for the processes along the dotted curve b (see Fig.l) are presented. 
The curve b corresponds to the case H — v + 1/a — 0.8. In the figure the y- axis indicates 
the values of s, whereas along the x-axis both the values of a (in the bottom, linear scale) 
and of v (in the top, non - linear scale) are shown. Black squares indicate the values of 
s measured along the curve b, whereas solid line indicates the value s = 0.8, which we 
expect to get. We see, that numerical results are well fitted by the expected line. For the 
comparison we show by crosses the exponent s measured from Eq. (18) for the second order 
structure function, q = 2. The theoretical second order structure function is infinite for the 
Levy motion. This circumstance has two consequences for the numerical simulation and/ or 
at experimental data processing. The first one is that the second order structure function 
increases with the length of the Levy motion trajectory increasing. The second consequence 
is that the exponent s does not depend on a and demonstrates "pseudo - Gaussian" behavior. 
Indeed, the dotted line s = v + 1/2 (which indicates the relation between s and v for the 
approximation to the fractional Brownian motion) well fits the relation indicated by crosses 
in Fig. 4b. 

In Fig. 5 we illustrate the studies of the range. We investigate r - dependence of the 
range R u of our approximation to the fLm with the Levy index a > 1, 
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R u (t)= sup [L a ,„{s) - L a>v {0)\ - inf [L a , v (s) - L a , v (0)] , (19) 

0<S<T 0<S<T 

(Wkt 5 ' . (20) 

Since the "t h law" is fulfilled for the range of the Levy motion, see Eqs.(6,7), it is expected, 
according our way of modeling, that H u is close to H — v + 1/a. 

In the empirical rescaled range analysis, that is, at experimental data processing, or in 
numerical simulation, the range of the random process is divided by the standard deviation 
of its increments after subtraction of a linear trend |pfl|| . This procedure, called the Hurst 
method, or the method of normalized range, in particular, smooths the variations of the range 
on different segments of time series. As the result of the empirical rescaled range analysis of 
experimental data, one gets the Hurst exponent of the process, that is, the exponent H v in 
our case. However, the Hurst method is not satisfactory for the Levy motion because of the 
infinity of the theoretical value of the standard deviation. Therefore, we propose to modify 
the Hurst method by exploiting the 1/a - th root of the a - th moment instead of standard 
deviation, that is, 

/ 1 T \iA» 

°« = [-Em)\ a ) (2i) 

Since it has only weak logarithmic divergence with the number of terms in the sum increasing, 
then the power - law dependence (20) is not changed. 

In Fig. 5a we present the results of applying the modified Hurst method to the ranges 
of the processes with the parameters along the vertical line a in Fig. 1. The Levy index 
a is equal 1.7. In the figure the exponent H u versus v is depicted by black squares. The 
relation H v — v + 1/a is shown by solid line. As in case of the structure function, see 
Fig. 4a, we notice that numerical results are well fitted by the line showing the expected 
relation. However, the discrepancy appear for v reaching its upper boundary, that is, for the 
strongly persistent case. In this respect we may repeat our above comment to the analogous 
discrepancies in case of the structure function. 

In Fig. 5b we present the results of applying the modified Hurst method to the ranges of 
the process with the parameters changing along the curve b in Fig. 1. This curve corresponds 
to the case H = u + l/a = 0.8. The y - axis indicates the values of H v , whereas along the x - 
axis both the values of a (in the bottom, linear scale) and of v (in the top, non - linear scale) 
are shown. Black squares indicate the values of H v measured along the curve b, whereas solid 
line indicates the value H v = 0.8, which we expect to get. We see the discrepancy between 
the numerical results and the expected line. However, the discrepancy is even larger, if one 
uses the ordinary, or non - modified, Hurst method. The results of its exploiting are shown 
by crosses, whereas the "pseudo - Gaussian" relation H v = u + 1/2 (which does not depend 
on a) is depicted by dotted line. One can clearly see the inapplicability of the "traditional" 
Hurst method for characterizing properties of the Levy motions. 

V. RESULTS 

The results of the paper are as follows. 
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1. We propose a model for the random process, whose increments are stationary, self 

- affine and distributed with the stable probability law. By analogy with the fractional 
Brownian motion, the family of these processes can be called fractional Levy motion. 

2. When constructing our model, the two basic steps are: 

(i) the use of the Gnedenko limit theorem for the normal attraction basin of a stable 
law; the theorem gives us a simple way for generating a sequence of independent stably 
distributed variables, which approximate a discrete - time white Levy noise; 

(ii) fractional integration/differentiation of a white noise; this procedure converts white 
noise into a fractional one, thus allowing us to approximate a discrete - time fractional Levy 
noise. 

3. We find the ranges of the change of the order of fractional integration/differentiation, 
inside of which the increments of the Levy motion approximation obtained possess the 
property of self - affinity This property manifests itself in so - called " t h laws" , that is, in 
the power - law time - dependence of both the structure function of the order q less than 
the Levy index a and the range. The relation between the exponent H and the order of 
fractional integration/differentiation is also obtained. 

4. We study the "t h law" for the structure functions of the approximation and find 
a good agreement between numerical results and the theory for the q - th order structure 
function, q < a. We also demonstrate that the second order structure function has a "pseudo 

- Gaussian" time behavior, which is irrespective of the Levy index a. 

5. We study the "r H law" for the range and demonstrate, that the "traditional" Hurst 
method exploiting the normalized rescaled range leads to a "pseudo - Gaussian" time - 
dependence, which is irrespective of the Levy index a. Thus, we propose the modified Hurst 
method, in which the 1/a - th root of the a - th moment of the sequence of increments is 
used instead of standard deviation. The modified Hurst method leads to the rescaled range 
time - dependence, which is much closer to the theory than that of normalized range. 

6. We conclude that our approximation is suitable for simulating persistent Levy motion 
with the Levy indexes varying between 1 and 2. We also discuss the possible reasons why 
our model works worse in the anti - persistent region and in the region of Levy indexes less 
than 1. 

7. As it concerns with the "pseudo - Gaussian" effects described above, they allow us 
to suggest that at estimating the second order structure function and normalized span from 
experimental data the "Levy nature " of them can be easily masked. This, in turn, rises 
an interesting task of developing statistical methods for extracting reliable characteristics 
from experimental data, for which Levy statistics can be expected from, e.g., some physical 
reasons. 
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Fig. 1. Shaded region indicates the admissible range of the order z/of fractional in- 
tegration/differentiation of a white Levy noise. The Levy index a is varied from 1 to 2. 
The vertical thick line at the right corresponds to the fractional Brownian motion, a — 2, 
— 1/2 < v < 1/2. The horizontal thick line v — corresponds to the ordinary Levy motion. 
Above the horizontal line the process is persistent, whereas below it is anti - persistent. 
Along the vertical dotted line a the Levy index a is 1.7. Along the dotted line b the relation 
v + 1/a — 0.8 holds. The self - similarity properties of the approximation to fLm along the 
dotted lines are illustrated in Figs. 4, 5. 

Fig. 2. Probability densities (above) and their asymptotics (below) are shown for the 
approximation to discrete - time white Levy noises with the Levy indexes a — 1.2 and 
a = 1.7. The probability densities and the asymptotics of the stable laws are indicated by 
solid lines. Below for the qualitative comparison the asymptotics of the Gaussian distribution 
are shown. 

Fig. 3. Above: typical samples of approximation to a discrete - time white Levy noise, 
v — (at the left) and fractional Levy noise, v = 0.3 (at the right). Below: corresponding 
trajectories of approximation to the ordinary Levy motion (at the left) and the fractional 
Levy motion (at the right). 

Fig. 4. (a) Exponent s of the structure function of order q — 1/4 versus the order of 
fractional integration u, see Eq.(18) (black points). Processes along the vertical dotted line 
a from Fig.l are investigated. The Levy index a is 1.7. Solid line indicates the relation 
s = v + 1/a. (b) Exponent s is indicated for the processes along the dotted line b from 
Fig.l (black points). The parameters a , v of these processes obey the law v + 1/a — 0.8. 
The order q of the structure function is 1/4. Solid line indicates the relation s = 0.8. The 
exponent s for the second order structure function is indicated by crosses. Dotted line shows 
the "pseudo - Gaussian" relation s — u + 1/2. 

Fig. 5. (a) Exponent H u for the range, see Eqs. (19), (20), versus v is depicted by 
black points. The numerical results are obtained with the use of the modified Hurst method 
proposed in the paper. Processes along the vertical dotted line a from Fig. 1 are investigated. 
The Levy index a is 1.7. Solid line indicates the relation H v — v + l/a. (b) Exponent H v is 
indicated for the processes along the dotted line b from Fig.l (black points). The parameters 
of these processes obey the law v + l/a = 0.8. The modified Hurst method is used. Solid line 
indicates the relation H u = 0.8. The exponent H v , obtained with the use of "traditional" 
Hurst method, is indicated by crosses. Dotted line shows the "pseudo - Gaussian" relation 
s = i/ + 1/2. 
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